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ABSTRACT 

We have investigated the brown dwarf (BD) and stellar disc fractions in the Upper 
Scorpius OB Association (USco) and compared them with several other young regions. 
We have compiled the most complete sample of of all spectroscopically confirmed BDs 
in USco, and have made use of the Wide Field Survey Explorer (WISE) catalog to 
identify the disc candidates. We report on the discovery of 12 new BD discs in USco, 
with SpT between M6 and M8.5. The WISE colors for the new discs are similar to 
the primordial (transition) discs earlier detected in USco. Combining with previous 
surveys, we find the lowest inner disc fractions ('^20-25%) for a wide range in stellar 
masses ('^O.O1-4.OM0) in the USco association. The low disc fractions for high-mass 
stars in USco (and the other clusters) are consistent with an evolutionary decline in 
inner disc frequency with age. However, BD disc fractions are higher than those for the 
stars in 1-3 Myr clusters, but very low in the ~5 Myr old USco. Also, primordial BD 
discs are still visible in the '^10 Myr old TW Hydrae association, whereas the higher 
mass stars have all transitioned to the debris stage by this age. The disc frequencies for 
BDs do not show any dependence on the stellar density or the BD/star number ratio in 
a cluster. We also find no convincing evidence that any of the well-known disc dispersal 
mechanisms for stars are active in BD discs. We suggest that the large differences in the 
observed BD disc fractions between regions may well be due to different BD formation 
mechanisms and therefore different initial disc fractions/properties. 

We also present a WISE SED classification scheme, based on the Kg and WISE 
bands of 3.4-12|j.m. We have determined certain thresholds in the WISE spectral slope 
versus spectral type diagrams to distinguish between the red population of Class 
I/II systems and the Class III sequence. We have found the WISE [3.4]-[12] color 
to provide the best distinction between the photospheric and the disc population. Our 
work includes a comparison of the sensitivities of WISE and Spitzer disc surveys. We 
estimate that WISE can be incomplete for discs at spectral type later than M8 in 
distant clusters such as SOri. WISE should be able to recover the M8-M9 discs in the 
nearby young clusters. 

Key words: circumstellar matter - stars: low-mass, brown dwarfs - open clusters 
and associations: individuals (Upper Scorpius OB Association) 



1 INTRODUCTION 

Recent surveys of nearby young clusters spanning a range 
in ages and masses have revealed extensive information on 
the disc frequencies and the disc dissipation timescales. The 
general picture of inner disc evolution that emerges from 
Spitzer/IRAC observations of several clusters shows a sharp 



decline in the disc frequency of young stars with ages be- 
tween ~1 and 3 Myr. By an age of ~5 Myr, the disc frac- 
tion has dropped by more than 80%, while the timescale for 
nearly all stars to lose their inner discs appears to be ~10 
Myr (e.g., Haisch et al. 2001; Hillenbrand 2005; Hernandez et 
al. 2007). Other than the overall disc fractions, the strength 
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in the disc emission at mid-infrared wavelengths also shows a 
decrease with age. A majority of young disc sources in Tau- 
rus (~1 Myr) possess optically thick inner discs, while most 
stars by an age of ~5 Myr possess "anemic" discs which dis- 
play weaker IRAC excesses than optically thick disc systems 
(e.g., Lada et al. 2006; Carpenter et al. 2009). The age of 
~5 Myr is thus critical, in the sense that nearly all optically 
thick primordial disc material in the inner ~1 AU region 
has dissipated by this age, and most discs begin to show a 
transition from a primordial to a debris disc system. 

The disc dissipation timescales, however, show a de- 
pendence on the stellar mass and the stellar density in the 
region, such that longer disc lifetimes are observed for the 
sub-stellar sources, and disc fractions are found to be higher 
in less dense environments (e.g., Riaz & Gizis 2008; Riaz et 
al. 2009). In Riaz et al. (2009), we had determined a BD disc 
fraction of ~10% in the ~5 Myr old Upper Scorpius OB asso- 
ciation (hereafter, USco). In comparison with younger ~l-3 
Myr old clusters that show a ~40-50% disc fraction, this in- 
dicated significant disc dispersal for sub-stellar sources over 
a short ~2 Myr time span. In contrast, the BD disc frac- 
tion in the ~10 Myr old TW Hydrae Association (hereafter, 
TWA; Webb et al. 1999) is much higher (50±20%; Riaz & 
Gizis 2008). We had argued that the BD density in the TWA 
is at least ~100 times less than in USco, which could prolong 
the disc lifetimes. Thus age may not be a dominant factor 
but a less dense environment could explain the higher disc 
fractions observed. 

In the present work, we have further investigated the 
dependence of the BD disc fractions on the age, density and 
the BD formation mechanisms in a cluster. We have focused 
again on the USco association, since it is at a critical inter- 
mediate age of ~5 Myr which is important to probe inner 
disc dissipation timescales. The USco association is part of 
the Scorpius Centaurus complex (de Geus et al. 1989), lo- 
cated at a distance of 145±2 pc (de Bruijne et al. 1997). 
The age of USco is ~5 Myr, with little scatter of ±2 Myr 
(e.g., Preibisch & Zinnecker 1999). The region is relatively 
free of extinction (Av ^ 2 mag) and star formation has al- 
ready ended (Walter et al. 1994). The association has been 
targeted at multiple wavelengths over the past decade start- 
ing with X-rays (e.g. Preibisch et al. 1998), with Hipparcos 
(de Zeeuw et al. 1999), in the optical (Preibisch et al. 2001; 
2002; Ardila et al. 2000; Martin et al. 2004; Slesnick et al. 
2006), and in the near-infrared (e.g., Lodieu et al. 2007), re- 
sulting in several low-mass stars and BDs being confirmed as 
genuine spectroscopic members (Ardila et al. 2000; Martin 
et al. 2004; Slesnick et al. 2006; Lodieu et al. 2008; Lodieu et 
al. 2011). The disc population in USco has been studied by 
Carpenter et al. (2006; 2009), Scholz et al. (2007) and Riaz 
et al. (2009). The B/A and F/G type stars are aU found 
to be debris disc systems, while primordial discs with inner 
disc clearing (no excess shortward of ~8vim) are more abun- 
dant among the K/M types and the sub-stellar sources. We 
have compiled a sample of all spectroscopically confirmed 
BD members in USco, and have made use of the recently 
released Wide Field Survey Explorer (WISE; Wright et al. 
2010) catalog to identify the disc candidates. Our sample 
also includes some M dwarfs in USco that have not been 
previously studied in the mid-infrared. In Section ^ we 
describe our WISE spectral energy distribution (SED) clas- 
sification scheme and the WISE colors that can best distin- 



guish between the discless and the disc population. Section 
33] presents the results on new disc sources found in USco, 
and Section t )4.2l compares the Ha emission strength with 
the disc emission. Section ^ discusses the various factors 
that may affect the disc fractions in a cluster. 



2 SAMPLE SELECTION 

We have included in our sample the optical [RI) survey 
complemented by 2MASS infrared photometry conducted by 
Slesnick et al. (2006; 2008) over the full association down to 
i? ~ 20 mag (corresponding to masses of ~0.02 M©). This is 
the most complete survey of the entire USco region. A large 
number of photometric candidates were selected and a total 
of 107 were confirmed as spectroscopic members based on 
strong Hq emission and weak gravity-sensitive features. The 
spectral types (SpT) for these sources are between M4 and 
M8.5. We note that the disc properties of 28 USco members 
with SpT later than M6 identified by Slesnick et al. (2006) 
are already discussed in Riaz et al. (2009) using Spitzer 
archival data. The USco association has been independently 
targeted by the UKIRT Infrared Deep Sky Survey (UKIDSS) 
Galactic Clusters Survey (GCS) in two optical filters (ZY) 
and three infrared passbands (JHK) (Dye et al. 2006). The 
depth of GCS in terms of mass is roughly 10 Jupiter masses. 
Lodieu et al. (2007) identified over 140 photometric candi- 
dates in 6.5 square degrees imaged by the GCS based on 
their photometry, colours, and proper motions. Over 90% of 
these candidates were confirmed spectroscopically as USco 
members using cross-dispersed near-infrared spectroscopy 
(Gemini/GNIRS; Lodieu et al. 2008) and muhi-fibre opti- 
cal spectroscopy (AAT/AAOmega; Lodieu et al. 2011). The 
SpT for these members range between M3.5 and M9. The 
GCS USco surveys are complete down to the SpT of M9. 

We have cross-matched this list of USco spectroscopic 
members with the WISE catalog. We adopted a matching 
radius of 6.5 arcsec to take into account the spatial resolu- 
tion of WISE, which is ^^6 times poorer than the UKIDSS 
seeing-limited images. We note that all WISE counterparts 
to the UKIDSS sources are within one arcsec (or less) as 
the result of the good astrometric precision of WISE (0.15 
arcsec; Wright et al. 2010). We will refer to the WISE 3.4, 
4.6, 12 and 22(j.m bands as WISE band 1, 2, 3 and 4, respec- 
tively. Our final USco- WISE matched catalog consists of 215 
sources and covers spectral types (SpT) between M3.5 and 
M8.5. We have considered the sources in our sample with 
SpT ^M6 to be BDs or sub-stellar sources, as described in 
Section ^X] Out of the 215 sources in our sample, 43 are BDs 
and 146 are M dwarfs. The USco samples previously studied 
by Carpenter et al. (2006; 2009) covered SpT of B0-M5. Our 
sample is complete for M4 and later types in USco. 

As discussed in Section ij3] we have considered the 
WISE 3.4, 4.6, and 12ji.m bands for SED classification, and 
have therefore required that all sources are detected at a 
>3-(T level in these bands. A majority of the USco sources 
lack detections in the WISE band 4 at a >3-a level. Our 
>3-(j quality criteria is to ensure good S/N detections, and 
to avoid spurious sources where the large excesses could be 
due to the background noise. This quality criteria imposed 
could introduce a potential bias to the disc fraction obtained, 
since the full sample detected in WISE at a >3-(j level would 
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be smaller in size than the number detected in e.g., a deep 
Spitzer survey. A disc fraction obtained with WISE could 
thus be higher than the "true" fraction inferred from the 
higher signal-to-noise Spitzer observations. This difference, 
however, is found to be small (~3%), and the disc fractions 
obtained with WISE in USco and the other comparison re- 
gions are found to be consistent within the uncertainties, as 
discussed in detail in Sections i|4.1l and 

What is the smallest excess emission that we could de- 
tect for the brown dwarfs with WISE? This can be deter- 
mined by calculating the ratio of the observed flux to the 
predicted photospheric flux, (Fobs/Fphot). As we discuss in 
Section the 12jxm point is vital in the confirmation 

of excess sources. The predicted photospheric flux at 12jxm 
and 22(i-m for an M6 dwarf is ^^0.25 mjy and ~0.I2 mjy, 
respectively. For the latest type objects in our sample at 
spectral type of M8-M8.5, the predicted photospheric flux at 
12(j.m and 22)i.m is ~0.15 mJy and ~0.04 mJy, respectively. 
The WISE 5-a point source sensitivity at 12nm and 22(j.m 
is ~0.7 mJy and ~5 mJy, respectively (WISE Explanatory 
Supplement). It is therefore impossible to detect the photo- 
spheres for the BDs at these wavelengths, and our survey in 
these bands is only sensitive to disc detections. If we con- 
sider the WISE 5-cr detection limit of ~0.7 mJy at 12nm 
to be the lowest detectable observed flux, then this implies 
Fobs/Fphot ~ 2.8 for the M6 sources and Fobs/Fphot ~ 
4.6 for the later type M8-M8.5 targets. Therefore, the mini- 
mum 12^.m excesses that we can detect will be at the level of 
Fobs/Fphot > 3-5 for our M6-M8.5 targets. At 22\xm, the 
minimum detectable excess would be Fobs/Fphot > 40 for 
the M6 sources and Fobs/Fphot > 125 for the later type tar- 
gets. These large limits would explain why most of our disc 
sources were not detected in the 22)i.m band (Section i |4.1|l . 
At 3.4[im and 4.6n.m, the predicted photospheric fluxes are 
between ~3 and 0.7 mJy for the M6-M8.5 spectral type 
range. These fluxes are higher than the WISE 5-a sensi- 
tivity of 0.04 mJy and 0.08 mJy in the 3.4)i.m and 4.6vi.m 
bands, respectively. Photospheric detection is thus possible 
in these shorter wavelength bands, and we should be able to 
detect any small excess above the photosphere. We note that 
the minimum excess that can be detected at any wavelength 
would also be dependent on the calibration uncertainty, the 
finite signal to noise, and the intrinsic spread in the stellar 
colors. To summarize, the two potential biases in our WISE 
disc survey are the quality criteria, and the minimum 

excess that can be detected at 12 and 22(im. We have dis- 
cussed the completeness level of our WISE disc survey in 
USco in Section >i4.1l and a comparison with the Spitzer 
surveys in other young regions is provided in Section 



3 WISE SED CLASSIFICATION 

In order to devise a SED classification scheme based on the 
observed fiuxes in the WISE bands, we started with the 
Spitzer/IRAC catalog of disc population in the Taurus star- 
forming region provided in Luhman et al. (2010). We cross- 
matched the Taurus sources in the WISE catalog and se- 
lected the ones that were detected at a level in WISE 
1, 2, and 3 bands. This Taurus- WISE catalog consists of 194 
sources. Luhman et al. (2010) have computed spectral slopes 
[a — d \og{\F\)/d log(A)] using four different wavelength 



ranges: (iS'./[8.0], is:./[24], [3.6]/[8.0], and [3.6]/[24]). They 
have then marked certain thresholds in the spectral slope 
versus SpT diagrams to distinguish between the red popu- 
lation of Class I and II sources and the sequence of Class III 
stars. For most of the sources, the 2-8(i.m {Ks/[8.0]) slope 
has been used for classification purposes. This particular 
spectral slope shows a rise in values towards later types of 
>M4. As mentioned, a majority of the USco sources lack 
detections in the WISE 4 band at 22fi.m. We have therefore 
computed the 2-12jxm (A's/[12]) spectral slopes to classify 
these sources. This wavelength range covers four data points 
at ~2.16, 3.4, 4.6 and 12)i.m. As we discuss later in this sec- 
tion, including the 12iJ.m point is essential as it is the only 
point that can identify the disc candidates. Nearly all of the 
USco sources show photospheric emission at shorter wave- 
lengths and most are undetected at 22jxm. We could use 
the 3.4-12fi.m wavelength range, although it would provide 
lesser number of points to determine the spectral slopes. 
We first computed WISE 2-12fi.m spectral slopes {awiss) 
for the Taurus sources, and compared these to the IRAC 
spectral slopes of 2-8jim from Luhman et al. (2010). Fig. [T] 
shows that the two slopes are well correlated, as indicated 
by the straight line fit to the plotted points. We have plot- 
ted with different colors the Taurus Class I/II/III sources. 
These are the SED classes assigned by Luhman et al. (2010) 
using the IRAC 2-8(xm slopes. The Class III sources in Fig.[T] 
have awiSE of 5C -2, the Class II sources are mainly con- 
centrated between awiSE of ~ -1.8 and ~ -0.2, while the 
Class I sources have awiSE ^ -0.2. The Class III sources 
show a more clear separation from the Class II systems. 
The boundary between Class I and II is less distinguished. 
Fig.[2]shows awiSE plotted against the SpT for the Taurus 
sources. There is a lack of sources with SpT later than M8 
that have good S/N detections in the 12(im band, and thus 
do not have awiSE estimates. The threshold of awiSE ~ 
-2 provides a good separation between the Class II and III 
sources. There may be a rise in this threshold at later types, 
such as the threshold may be at -1.5 at SpT of M8.5 and 
later. However, we do not have enough points at later types 
to confirm a rise in the Class II/III threshold. The threshold 
between Class I and II is around ~ -0.2, although a higher 
value of awiSE^ 0.0 is required to distinguish between these 
two classes at earlier types. 

In Fig. [S] we have compared the WISE spectral slopes 
with various WISE and WISE-NIR colors for the Taurus 
sources. For colors including the 22(i.m photometry, we have 
only plotted the objects with a S/N^3 in the 22(im band. 
The [3.4]-[12] and _K'i,-[12] colors {top panel) are well cor- 
related with the WISE slopes and can best distinguish be- 
tween the different SED classes. We can consider a boundary 
of [3.4]-[12] ~1.2 and K^-[12] ~1.8 to separate the Class III 
sources from the disc population of Class I and II systems. 
The color [4.6]-[12] {2nd panel) also provides a good sep- 
aration between Class II and Class III sources (boundary 
of [4.6]-[12]~l). The color [3.4]-[4.6] is vague in separating 
the Class II and III sources, but can be used to identify the 
protostellar systems ([3.4]-[4.6]>l). The bottom two pan- 
els in Fig. [3] plot the WISE colors which include the 22(i.m 
photometry. The typical range in these colors for the photo- 
spheric sources is between and 1. However, there are a few 
evolved/transition discs among the Class HI objects which 
show redder colors between '^l and 3. A transition disc ex- 
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hibits photospheric emission at wavelengths shortward of 
'^lOum due to inner disc clearing, but retains excess emis- 
sion at longer wavelengths due to the presence of circumstel- 
lar material in the outer disc regions. Such transition discs 
can be more clearly separated from the main Class III locus 
in the [4.6]- [22] and the [12]- [22] colors (boundary of >1). To 
summarize, colors which include the 12(i.m magnitudes can 
more clearly make a distinction between the photospheric 
and the disc population, whereas the colors which include 
the 22|i.m magnitudes can be important in identifying the 
transition disc candidates. 

Based on the analysis above, we have classified the USco 
sources as Class III if awiSE ^ -2.0, and as Class II if -2.0 
<CiwiSE< -0.2. Figure |4] shows the SpT vs. awiSE dia- 
gram for the full USco sample. The disc sources identified 
are above the threshold of awiSE = -2.0, and have been 
classified as Class II systems. The extinction estimate in 
USco is Ay < 2 mag (Walter et al. 1994). We have therefore 
not applied any reddening corrections to the spectral slopes. 
The values for otwiSE and the SED classes assigned to the 
USco sources are listed in Table [1] 



4 RESULTS 

4.1 New Discs in USco 

Out of a sample of 215 sources, we have found 48 new discs 
in USco, 12 of which are BD discs. The SEDs for the disc 
sources are shown in Figs. [7|and[51 To fit the stellar photo- 
sphere, we have considered the NextGen models (Hauschildt 
et al. 1999) for solar metallicity and log g = 3.5. The ef- 
fective temperature was obtained from the SpT using the 
relation for young sources in Luhman et al. (2003a). We 
have checked the various fiags associated with the WISE 
photometry. Out of the 48 disc sources, 40 have 'cc_fiag' of 
'0000', indicating that these are not affected by any artifacts. 
The other disc sources either have cc_flag=hOOO or dOOO, im- 
plying that the photometry in band 1 may be affected by 
the halo of a nearby bright source or a diffraction spike. 
However, we have checked the WISE and UKIDSS images 
and do not find the sources to be extended or blended with 
any bright nearby star, or a spike located close to or across 
the object. The photometric quality flag, ph_qual, for these 
sources is either 'AAAA' or 'AAAB'. Thus even if band 1 
may be affected by some artifact, the photometric quality 
is the 'highest', so these should be reliable detections. One 
source (UScol62847) is blended with a bright nearby star in 
bands 1, 2 and 3 (cc_fiag is 'h' in these 3 bands). This object 
has ph_qual = 'AAAB', indicating 'highest' quality photom- 
etry in bands 1, 2, and 3. This source is resolved into 2 bright 
stars in the UKIDSS GCS images but lies close to a nearby 
bright star in the WISE images. The WISE resolution is 6.5" 
which is large compared to seeing-limited images taken by 
UKIDSS. Since the composite source shows an excess emis- 
sion at 12 and 22nm, we have considered this source as a 
disc candidate. The WISE photometry for the disc sources 
is listed in Table[Tl The value of '99.999' in Table [T] indicates 
a non-detection in that band. 

In Fig. O we have compared the WISE colors for the 
disc sources from our work with the primordial and debris 
discs in USco found by Carpenter et al. (2009). The range in 



uncertainties are 0.03-0.06 mag for the ii's-[3.4] color, 0.03- 
0.15 mag for the [3.4]- [12] color, and 0.04-0.2 mag for the 
[4.6]-[22] color. The colors for all primordial discs are much 
redder than the debris disc sources. Carpenter et al. have 
noted a high fraction of primordial discs among the M-type 
stars, the K/G types are mainly all debris discs with a few 
primordial discs, while the higher mass B/A/F stars are all 
debris disc systems. Their categorization of a primordial disc 
is based on excess emission observed at 8 and 16[im, while a 
debris disc is classified based on excess emission at ^24|i.m 
(Section Thus these primordial discs can be considered 
as primordial transition discs. The age of ~5 Myr for USco is 
considered to be a transitional age, when circumstellar discs 
are in a transition from a primordial to a debris disc system. 
A majority of disc sources at this age are found to exhibit in- 
ner disc clearing. In the literature, different definitions have 
been used to distinguish between a primordial transition disc 
and a debris disc. A primordial transition disc is basically 
defined to be a disc that shows photospheric emission short- 
ward of ~10(i.m, indicating that the inner ~0.1-1 AU disc 
region (depending on the mass of the central source) is de- 
void of optically thick material, but it retains strong excess 
emission at longer wavelengths of >10vim, which indicates 
the presence of primordial material in the outer disc regions. 
The color excesses at longer wavelengths for the primordial 
transition discs would be comparable in strength to the main 
Class II population in the cluster. 

Fig. [5] shows that the WISE colors for our M dwarf 
and BD discs are similar to the primordial transition 
discs in USco. If we consider the Luhman et al. (2010) 
boundary of K-[3.6] ~0.4 for photospheric colors of M4- 
M5 stars, and K-[3.6] ~0.5 for M6-M8 objects, then we 
have about 11 M dwarf discs and 3 BD discs that lie 
above this boundary. Among these sources, we have two 
BD discs, UScol61005 and UScol62847, and 6 M dwarf 
discs, UScol55140, UScol61117, UScol61349, UScol62222, 
UScol61816, and UScol61819, which have Fobs/Fphot ~2- 
2.5 at 3.4 and 4.6(i.m, while the rest have Fobs/Fphot ~ 1.1- 
1.4 (Figs.[7l[8]|. A large fraction of USco discs thus show signs 
of inner disc clearing, with weak or no excess emission short- 
ward of 12)i.m. We have therefore classified the new discs 
found in our survey as primordial transition discs. Among 
the BD discs, four sources, UScol55420 (M8), UScol60953 
(M6), UScol61940 (M6) and UScol63027 (M6), show flaring 
between 12 and 22(i.m (Fig. [T]). For most of the previously 
known USco brown dwarf discs, the 24|i.m flux density is be- 
tween ~2 and 5 mjy, and only a few discs have been found 
with higher fluxes of 15-20 mJy at 24[im (Scholz et al. 2007, 
Riaz et al. 2009) . A ~2-5 mJy flux density will be below the 
WISE 5-a point-source sensitivity of 6 mJy at 22[j.m (WISE 
Explanatory Supplement), which could explain why most of 
these discs were not detected in the 22\xm band. In Fig. [B] 
the USco BD discs show larger color excesses compared to 
the earlier type discs. The 12(i.m to Ks flux ratio for most 
BD discs is also larger in comparison to earlier type discs, 
indicating larger 12[j.m flux excess for the BD discs. Since 
these disc sources are presumably at the same age, this sug- 
gests that discs persist for a longer time around later type 
stars. We discuss this further in Section ^ 

All of the disc sources show excess emission at 12[i.m, 
including the BD discs, and there are no discs that have ex- 
cesses at 22[j.m but none at 12(i-m. Based on WISE data, a 
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debris disc can be defined as a source with excess at 22(j.m 
only. The minimum excess that we can detect with WISE for 
a brown dwarf disc at 22(i.m would have to be at least ~40 
times larger than the photospheric flux at this wavelength. 
Therefore WISE may have missed brown dwarf debris discs 
due to low sensitivities. However, we note that none of the 
deeper Spitzer surveys have found BD debris discs, i.e. discs 
that showed excess emission at 24nm but none at shorter 
wavelengths. In comparison, for an M4 dwarf, an excess ~5 
times larger than the photospheric flux could have been de- 
tected with WISE. As mentioned, a high fraction of M dwarf 
discs in USco are found to be in the primordial stage, com- 
pared to earlier-type stars (Carpenter et al. 2009). The ab- 
sence of debris systems among our M dwarf and BD targets 
is thus consistent with previous results. 

From our present survey in USco, we flnd a BD disc frac- 
tion of 28±9% (12/43). The error bars are the 1-a Gaussian 
distribution uncertainties. This fraction is comparable to the 
37±9% fraction reported by Scholz et al. (2007) from their 
~8-14[im Spitzer/IRS spectral observations. In Riaz et al. 
(2009) , we had reported a much lower USco BD disc fraction 
of 11±4% (3/28). This fraction was based on Spitzer /IRAC 
and MIPS 24)i.m observations, and was mainly determined 
from the excess emission observed at 8[i.m. Nearly 90% of the 
sources in that survey were undetected at 24 nm. As we have 
found in the present work, a majority of the discs show pho- 
tospheric emission shortward of 12(i.m. A small disc fraction 
based on S^m excess emission is thus expected. Combin- 
ing our present work with Scholz et al. (2007) and Riaz et 
al. (2009), we have a USco BD disc fraction of 26±9%. For 
the M dwarfs in USco, we have a disc fraction of 25±8% 
(36/146). Our M dwarf disc fraction is comparable (within 
the uncertainties) to the 17±6% fraction reported by Car- 
penter et al. (2009) for the M type primordial discs in their 
sample. 

We note that due to lower sensitivities, a disc survey 
conducted with WISE could be incomplete as compared to 
earlier deep Spitzer surveys. There are 16 known brown 
dwarf discs in USco detected by previous Spitzer surveys 
(Scholz et al. 2007; Riaz et al. 2009). The 1-a sensitivity at 
8-lliJ.m from these Spitzer surveys is ~0.09 mjy. These sur- 
veys were thus deep enough to detect the photospheres for 
the brown dwarf targets. Out of these 16 Spitzer BD discs, 
15 were detected in the WISE 1, 2, and 3 bands at S/N > 
3. One of the discs (usdl61916) was detected in the WISE 

1 and 2 bands, but was undetected at 12)i.m. The observed 
flux density for this object at lljim is 0.65 mJy (Scholz et 
al. 2007), which is below the WISE 5-a sensitivity at 12nm 
(^0.7 mJy). This is an M8 disc. However, the M9 discs from 
Scholz et al. (2007) Spitzer survey were also recovered by 
WISE at a >3-(j level in the 1, 2, and 3 bands. Half of 
these Spitzer discs are detected in the WISE 22(i.m band at 
a S/N>3, while the rest are undetected at this wavelength 
(S/N < 2). Other than the WISE lower sensitivities, another 
potential bias could be due to the quality criteria that we 
have applied, in terms of the S/N to be at least 3 in bands 1, 

2 and 3, which would reduce the size of the original sample 
and bias the results towards a higher disc fraction. We have 
checked for WISE matches for the full Spitzer sample of 63 
sources, including the discless objects, from the work of Riaz 
et al. (2009) and Scholz et al. (2007). There are 11 objects 
which have S/N < 3 in the WISE 12nm band, while the rest 



were detected at higher S/N in all three WISE bands. The 
BD disc fraction in USco from Spitzer surveys was 16/63 = 
25.4%± 9%. Based on the WISE-recovered sample, the BD 
disc fraction wiU be 15/52 = 28.8% ± 9%. These fractions 
are still consistent within the uncertainties and the differ- 
ences are not significant. We can therefore estimate that the 
lower WISE sensitivities and our quality criteria would re- 
sult in a BD disc fraction in USco which is slightly higher 
(by ~3%) than the fraction determined from the higher S /N 
Spitzer data. 



4.2 Hq versus Disc Emission in USco 

Figure [9] compares the Ho emission strength for the disc 
sources in USco. The new discs we have found in this work 
are marked by black squares. The Ha equivalent width and 
SpT are from Slesnick et al. (2006; 2008) and Lodieu et 
al. (2008; 2011). We have included in this figure previously 
known USco discs from Scholz et al. (2007) and Riaz et 
al. (2009). The dashed line in this figure is the empirical 
accretor/non-accretor boundary from Barrado y Navascues 
& Martm (2003). Sources that lie above this boundary can 
be considered as CTTS, and the ones below it as the non- 
accreting weak-line T Tauri stars (WTTS). We find a nearly 
equal number of disc sources above and below the accretor 
boundary, which indicates that the presence of passive or 
non-accreting discs is as common in USco as accreting disc 
sources. The disc bearing sources among the BDs (SpT^M6) 
are mostly all accreting systems. Emission in Hq is known 
to rise towards the later types, with a peak near SpT of 
M6-M7 (e.g., Riaz et al. 2006). It is thus more likely to find 
disc sources among the later types that also show strong 
emission in Ha. 

A lack of correlation between Class II and CTTS sys- 
tems, or Class III and WTTS systems, is expected at this 
transition age of ~5 Myr. A large fraction of USco discs 
show photospheric emission shortward of 12|i.m, indicating 
the inner ~0.1-1 AU disc region to be cleared of circum- 
stellar material. The Ha emission line traces the gas in the 
inner disc regions and is not a diagnostic of the total gas 
content in the disc. It thus seems reasonable to find an anti- 
correlation between the presence of disc emission and lack 
of Ha emission, since the disc emission is being detected 
at longer wavelengths and not from the inner part of the 
disc close to the central source. Fig. [9] thus indicates that 
selecting young stellar objects based on the strength in Ha 
emission could result in a rejection of young disc sources 
with inner disc clearings. 



5 DISCUSSION 

5.1 Disc Fractions versus Stellar Mass, Age and 
Density 

In Fig. 1101 we have plotted the disc fractions for clusters at 
ages of ~1-10 Myr in 3 mass bins: the sub-stellar sources 
(SpT=M6-M9; <0.1 Mq), the low-mass stars (SpT=K5- 
M5; ~0.1-1 Mq), and the high-mass stars (SpT<K5; ~1- 
4 Mq). Sections]^ provides details on the compilation of 
the disc fractions as well as the fraction of discs that WISE 
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might have missed around M6-M8.5 type objects due to lim- 
ited sensitivity. The disc fractions plotted here imply the 'in- 
ner' disc fractions, since these have been obtained from ob- 
servations at ~3-24nm, which probe the inner disc regions of 
radii within ~1 AU in BD discs and within ~5 AU in higher 
mass T Tauri discs (e.g., Hillenbrand 2005). It is impor- 
tant to determine the evolutionary stage of these discs, i.e. 
whether these are in the primordial, primordial transition, 
or debris phase. As discussed in detail in Section ^X] nearly 
all of the discs in Taurus, Cha I, IC 348 and SOri are in the 
primordial/primordial transition stage, and a small ~l-2% 
fraction is found to be of debris discs. On the other hand, a 
majority of the high- and low-mass discs in USco and TWA 
are in the debris phase, while the BD discs in these older 
clusters can all be classified as primordial/primordial tran- 
sition discs. Thus the later type discs tend to remain in the 
primordial stage for a longer time. We note that the uncer- 
tainties on the disc fractions are quite large for some of these 
young clusters, which would make the fractions comparable 
to each other. Nevertheless, we find some obvious differences 
in the disc frequency trends for the different mass bins, as 
discussed further below. In Fig. 1111 we have plotted the disc 
fractions versus the stellar density in a cluster. By stellar 
density, we imply the total number of objects (stars -I- BDs) 
detected in a given survey area. Details on the compilation 
of stellar densities are provided in SectionfjX] Our results for 
stellar densities indicate SOri and Cha I to be denser than 
Taurus and USco, which is in agreement with what has been 
previously noted by several surveys. 

The disc fractions in Fig. [TO] are the lowest and nearly 
constant in USco over a wide range in stellar masses, unlike 
the SOri or Taurus clusters. The disc fractions decrease with 
stellar mass in Taurus, while the opposite trend is observed 
in SOri. IC 348 also shows a slight increase in disc fractions 
with decreasing stellar mass. In comparison, Cha I which 
is at a similar age of ~2 Myr as IC 348, shows a higher 
disc fraction for the high-mass stars. Luhman et al. (2010) 
had noted that IC 348 and SOri are denser than Cha I or 
Taurus, and that a denser environment could result in a 
lower disc fraction for the high-mass stars in a cluster. This 
can be seen in Fig. [11] (top panel), where the fractions for 
the high-mass stars in the denser regions of SOri and IC 
348 are lower than in Taurus or Cha I. However, the disc 
fractions are lower in USco than SOri, even though the two 
clusters are at similar ages (~3-5 Myr) and USco is even 
more sparse, with a stellar density ~20 times lower than 
in SOri. A denser environment also does not seem to affect 
the low-mass and the sub-stellar sources (Fig. Ilip . The disc 
fractions for these mass groups are similar in IC 348 and 
Cha I, which have about an order of a magnitude difference 
in stellar densities. Pfalzner et al. (2006) have simulated the 
relative disc mass loss induced by encounters between stars 
in a cluster. At an age of 2 Myr, their simulations indicate a 
disc mass loss of ~30% for intermediate-mass (1-10Mq) and 
low-mass (<1Mq) stars, whereas for stars with M >10Mq, 
it is --^60% or higher. As they explain, the average number 
of stellar collisions for high mass stars (>10Mq) is higher 
compared to intermediate- or low-mass stars, and so the 
disc dispersion rate will also be higher for the high mass 
sources. Also, if the disc sizes for BDs are smaller than high- 
mass stars, then frequent fiyby events would be less likely 
to disrupt the sub-stellar discs. While the mass ranges from 



Pfalzner et al. are different than the ones considered here, 
their simulations could explain why we see similar low-mass 
and BD disc fractions in Cha I and IC 348, compared to 
a factor of ~3 difference observed for the high-mass stars 
fFig. Ilip . A denser environment thus may be more effective 
in increasing the disc disruption rate for the high-mass stars. 

An alternate explanation for the prominent difference 
between the Taurus and SOri/USco disc fractions for the 
high-mass stars could be their early transition from a pri- 
mordial to a debris phase. By an age of ~5 Myr, ~86% of the 
high-mass discs are in the debris phase. Thus nearly all mas- 
sive sources have experienced significant inner disc clearing 
by this age, and are found to have excesses only at wave- 
lengths of ^24n.m, with photospheric emission at shorter 
wavelengths. The disc fraction then shows a rise again at 
older ages, from ~20% at 5 Myr to ~50% at 10 Myr. All 
of these older discs are in the debris phase, as these are 
mainly detected due to strong 24^-m or 70n.m excess emis- 
sion. Debris disc particles undergo continual shattering col- 
lisions with planetesimals and larger bodies in the system, 
due to which debris dust is continually generated. This re- 
sults in an increase in the dust mass at the onset of the 
debris disc phase. The increase in the debris disc fraction 
is due to the detection of large mass of cold dust in the 
outer regions of the disc, while the inner warm regions have 
been dissipated. We note that the TWA disc fraction has a 
high uncertainty due to the sparse population of this region 
(Section However, Siegler et al. (2007) and Reike et al. 
(2005) have reported disc fractions of ~40-50% for B-A type 
stars in 30-50 Myr old clusters, based on much larger sam- 
ple sizes. To summarize, the dip in the USco fraction for the 
high-mass stars is due to the decline of primordial discs by 
this age, or the fewer number of discs detected at this age 
with excesses at short wavelengths (s;24jim). 

In comparison to the high-mass stars, BD discs seem to 
evolve at a slower rate, as their inner disc fraction remains 
nearly constant until ~3 Myr. Thereafter, the disc fractions 
show a steep decline by a factor of ~2 between 3 and 5 
Myr. The statistical significance of this sharp drop in frac- 
tions between 3 and 5 Myr is discussed in Section [jX] This 
rapid inner disc dispersal over a short ~2 Myr timescale thus 
seems to occur for both the high-mass and the sub-stellar 
sources, albeit at an earlier age of 2-3 Myr for the high-mass 
stars, but at an older age of 3-5 Myr for the BDs. The later 
drop for the BDs could be explained by their slower viscous 
evolution timescales. The viscous scaling time for BDs are 
expected to be of the order of 10^ yr, longer than the typi- 
cal 10^-10^ yr for T Tauri stars (e.g., Alexander & Armitage 
2009). In other words, if accretion of material onto the BD 
occurs at a slower rate than higher mass stars, then discs 
could persist for a longer timescale. The inner discs around 
higher mass stars may also be replenished faster by material 
coming in from the outer disc region. Such a process of outer 
disc accretion seems to occur either at a slower rate for BDs, 
or may not be that significant if the disc sizes are too small 
compared to the higher mass T Tauri discs. It is interesting 
to note that the high mass and sub-stellar sources provide 
two extreme cases of rapid inner disc dispersion, while the 
intermediate low-mass stars show a more linear decline in 
disc fractions with age (Fig. llOp . These sources have an al- 
most 50/50 debris/primordial disc fraction at 5 Myr, which 
increases to 70/30 by 10 Myr (Table [2]). The low- mass discs 
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also show a rise at ages older than 10 Myr, due to an in- 
crease in the debris disc fraction (e.g., Siegler et al. 2007). 
Thus a large fraction of the low-mass discs are in the debris 
phase by ~10 Myr, but the transition is more gradual than 
that seen for the high-mass stars, with no sudden drops in 
the disc fractions observed. 

The rapid inner disc dispersal could be due to processes 
such as dust coagulation or planetesimal formation. The pro- 
cesses of significant grain growth and dust settling in the in- 
ner disc regions could lead to rapid formation of kilometer- 
sized large bodies or planatesimals in the disc. From the sim- 
ulations of DuUemond & Dominik (2004), such processes oc- 
cur on a faster timescale in the inner disc regions, and small 
grains can be removed from the inner regions on a timescale 
of less than 1 Myr. Photoevaporation could also result in 
rapid inner disc clearing. The process of photoevaporation 
sets in when the disc accretion rate falls below the pho- 
toevaporative wind rate, which is about W'^^Mq/jt (e.g., 
MuzeroUe et al. 2006). From the simulations of Alexander et 
al. (2009), giant planet formation and photoevaporation can 
both dominate inner disc clearing. At early times, planet for- 
mation/migration is the more dominating process as the ac- 
cretion rates are too high for photoevaporation to set in. At 
ages of ^2 Myr, the disc accretion rate falls below the pho- 
toevaporative wind rate and essentially all discs are being 
cleared out by photoevaporation. By an age of ~5 Myr, the 
inner disc fraction in their simulations has dropped to ~20% 
from nearly 100% at ~2 Myr. Thus a coupling of the planet 
formation/migration and photoevaporation processes could 
result in rapid inner disc clearing over a ~2-3 Myr timescale. 
Such processes would work for the high-mass stars where 
accretion rates measured from UV continuum emission are 
an age of ~1 Myr (e.g., Herczeg et 
al. 2008). For BDs, the accretion rates are < lO"^°M0/yr 
at an age of ~1 Myr (e.g., MuzeroUe et al. 2005), which 
would imply that photevaporation can be effective even at 
these early ages. The inner disc fractions for BDs at ~l-3 
Myr should then be as small as the fraction found in the 
USco. Giant planet formation is also feasible for large disc 
masses of ~0.1 Mq as considered by Alexander et al. (2009), 
whereas typical disc masses for BDs in Taurus are found to 
be much smaller, of the order of ~10~* Mq (e.g., Scholz et 
al. 2007). From the simulations by Payne & Lodato (2007), 
a close-in rocky planet with a mass of ~0.1 Mq can form 
around a BD at a semi-major axis as small as ~0.2 AU, 
considering an average BD disc mass of 1.5 xIQ-^Mq and 
an outer disc radius of 100 AU. However, the growth of such 
planetary cores is at a much slower rate around BDs com- 
pared to the solar-type stars. Dust coagulation leading to 
planet formation thus would take longer for the sub-stellar 
case, and cannot explain the sharp drop observed between 
the SOri and USco inner disc fractions in Fig. 1101 

Other processes such as photoionisation could also be 
in play. The ionizing radiation and winds from massive stars 
can affect their immediate surroundings, thereby dispersing 
the circumstellar material around stars. In Fig. 1121 we have 
plotted the known OB stars in USco, along with the disc 
sources from all mass groups. We do not find any particu- 
lar concentration of the disc sources farther away from the 
massive stars, which might have suggested minimal inter- 
action with the ionizing winds and so a higher probability 
of disc survival. Among the low-mass stars, there is a small 



concentration of disc sources near {a, (5)~(242°, -22°). Also, 
nearly all of the BD discs are located south of 5 ^ -21°, 
while most of the higher mass debris discs are located north 
of this declination. There may be a possible mass depen- 
dency in the spatial distribution of the discs. We note that 
these discs come from different surveys that may not cover 
the same region in USco. Overall, the discs are largely spread 
throughout the association, and we do not find any strong 
evidence of photoionisation in USco. Also, both SOri and 
USco are known to host OB stars. We would thus expect 
similar disc fractions in these regions if photoionization is 
the main cause of inner disc dissipation. 

To summarize the discussion above, we can explain the 
rapid decline in the high-mass inner disc fraction between 
~2 and 3 Myr using planet formation and/or photoevap- 
oration models. Accretion of material onto planets as well 
as photoevaporation would lead to a transition from a gas- 
rich optically thick disc to a secondary debris disc. This 
would explain the high fraction of debris disc systems found 
in USco for the high- mass stars (Table [2]). Stellar density 
may also play a role and disc lifetimes may be prolonged in 
a less dense environment. But this is only notable for the 
high-mass stars. None of the processes above can explain 
the drop observed for the BD disc fractions between 3 and 5 
Myr. In order to determine the significance of the observed 
decline between 3 and 5 Myr, we have built a probability 
distribution for the USco BD disc fraction using the method 
described in Burgasser et al. (2003). The distribution was 
constructed using a sample size of N=106, and number of 
discs n=28 (combining this work with Scholz et al. and Riaz 
et al. surveys). The peak is at 26.4% with 1-sigma Gaussian 
uncertainties of 9%. The 3 Myr SOri disc fraction is 57%± 
19%, i.e. the disc fraction could be between 38% and 76%. 
From the probability distribution, the probability that the 
USco sub-stellar disc fraction is between 38% and 76% is 
found to be ~28%, at a confidence level of 97.5%. There is 
thus a ~70% probability that the 3 and 5 Myr disc frac- 
tions are different. We note that the USco point is not just 
1-sigma below the SOri point, it is also 1-sigma lower than 
the Cha I, IC 348 and Taurus points. The question we have 
tried to address in this paper is not how likely it is that 
USco could actually be the same as any *one* of the others, 
but that it could be the same as *all* of the other clusters, 
which seems unlikely. Also notable for the BD discs is their 
persistent optically thick emission even at later ages. By ~5 
Myr, ~86% of the high-mass sources are in the debris phase, 
while for the BDs, even the ~10 Myr old discs show excesses 
at wavelengths as short as ~5|i.m, and are in their primor- 
dial phase (Riaz & Gizis 2008). In Fig. 1101 comparing the 
SOri and TWA disc fractions for BDs suggests that it may 
be nearly constant at ~50% over 1-10 Myr ages. The sharp 
drop between 3 and 5 Myr may not be a real trend with age, 
but perhaps the intrinsic or the original BD disc fraction in 
USco is much smaller than the younger or the older clusters. 

5.2 BD Formation Mechanisms 

Could it be possible that this particular low BD disc fraction 
in USco is due to the different formation mechanism of BDs 
dominant in these clusters? As noted earlier, Taurus, Cha 
I and USco are low density clusters and always have been, 
but IC 348 and SOri are more classic denser clusters. It may 
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be that the denser clusters formed more BDs from collaps- 
ing low-mass cores rather than other mechanisms such as, 
by fragmenting massive discs or ejection of protostellar em- 
bryos, and so initially they had a higher disc fraction. In 
Fig. 1131 we have plotted the BD disc fraction versus the ra- 
tio of the number of BDs to the number of stars in these 
clusters. Details on the compilation of the BD/Star number 
ratios are provided in Section fjX] USco has a factor of ~2 
higher ratio than the other clusters but the lowest disc frac- 
tion, which indicates that a higher number of BDs formed 
does not necessarily translate into a higher disc fraction. 
Also, the denser clusters such as SOri do not form a higher 
number of BDs, as the ratios are quite similar to the less 
dense Taurus and Cha I. 

Different formation mechanisms of BDs in these clusters 
could also lead to different disc properties, even if the total 
number of BDs formed is similar. There are five standard 
models considered for the formation of BDs: 

(i) From very low-mass cores, which should result in the 
presence of a large disc around the sub-stellar source formed 
(e.g., Padoan & Nordlund 2004; Hennebelle & Chabrier 
2008; Machida et al. 2009) 

(ii) Ejection from larger cores, which would result in a 
small, truncated disc at best (e.g., Reipurth & Clarke 2001; 
Goodwin et al. 2004) 

(iii) From massive discs; it is estimated that ~20-70% of 
the BDs formed by disc fragmentation can have discs (e.g., 
Stamatellos & Whitworth 2009ab) 

(iv) Liberated from binaries, which should result in the 
presence of a disc (e.g., Goodwin & Whitworth 2007) 

(v) In small 'cores' in filaments, the BD disc fraction 
should be the same as the low- mass disc fraction (e.g., Bon- 
nell et al. 2008; Clarke et al. 2008; Bate 2009) 

To clarify these models further, model (1) is based on a 
star-like formation of BDs, i.e. via the collapse of molec- 
ular cloud cores with sub-stellar masses. In this model, 
gravitationally unstable protostellar cores of BD mass are 
formed directly by the process of turbulent fragmentation 
(e.g., Padoan & Nordlund 2004). Since this is a star-like 
formation, and since most stars are formed with circumstel- 
lar discs, BDs formed via this mechanism are expected to 
be surrounded by discs. Model (2) considers that BDs are 
not formed through the collapse of low-mass cores, but that 
rather they are ejected as a result of encounters in multiple 
systems made up of a small number of stellar embryos. The 
ejected BDs have suffered multiple encounters that will have 
pruned their circumstellar discs. As a result, this model pre- 
dicts weaker accretion diagnostics and few BDs with discs 
(e.g., Reipurth & Clarke 2001). In model (3), BDs are formed 
due to gravitational fragmentation of massive extended discs 
around Sun-like primary stars. At an early age of ~20 kyr, 
this model predicts 23% of the BDs liberated in the field 
to have discs. Many of the BDs that still remain bound to 
the central star are loosely bound, and so when they escape, 
they are expected to retain their discs. The disc fraction 
could thus increase up to a maximum of 70% (Stamatel- 
los & Whitworth 2009ab). Model (4) is based on a binary 
disruption mechanism. In this model, BDs form as distant 
companions to low-mass stars, in particular M-dwarfs. Such 
systems are then readily disrupted by the mild perturba- 
tions of passing stars at relatively large distances (a few hun- 



dred AU or more). The result is a population of single BDs 
and low-mass hydrogen- burning stars. Since the disruption 
is gentle, both the BDs and the low-mass hydrogen-burning 
stars are able to retain their circumstellar discs (e.g., Good- 
win & Whitworth 2007). In model (5), BDs form in filamen- 
tary structures that develop as the gas falls in the gravita- 
tional potential of the forming cluster. In this model, there 
is no need for close interactions or ejections to ensure the 
formation of sub-stellar cores. Ejections can occur but are 
not fundamental to the process. This model predicts BDs to 
have the same disc properties as the low-mass stars formed, 
i.e., the circumstellar disc properties for BDs should form 
a continuum with low-mass stars. However, this mechanism 
would only be expected to be efficient in dense clusters. The 
abundance of BDs in stellar clusters is of the order of 25% 
in highly clustered regions, whereas it decreases to the order 
of 10% in isolated regions (Bonnell et al. 2008; Clarke et al. 
2008). 

Any, some, or all of these mechanisms could operate in a 
cluster, perhaps with the same relative efficiency, or perhaps 
different in different regions. Different mechanisms should 
give different disc properties relative to the other mecha- 
nisms, for e.g., very low-mass cores would probably produce 
fairly massive extended discs, while ejection would produce 
lower-mass truncated discs, or no discs at all. Once formed, 
we can probably assume that the same discs evolve in the 
same way independent of the formation mechanism. But the 
initial disc fraction may be different and discs could start sig- 
nificantly less massive/extended depending on the formation 
mechanism. The observations seem to suggest that BD discs 
do not have a nice decline in fraction with time. If different 
clusters form BDs via different mechanisms, then the initial 
BD disc populations will be different in these clusters. The 
observed BD disc fraction thus may only partially he due to 
age evolution, and rather the disc fraction may he a signature 
of the BD formation mechanism dominant in the cluster. For 
e.g., the mechanism (2) may be more active in USco than 
the rest. The ejection model (2) could result in an excess 
of BDs being formed, but very few of those BDs will have 
circumstellar discs. This mechanism may also be prevalent 
in Taurus that has a lower BD disc fraction than SOri de- 
spite being younger than this region. Model (3) may also be 
applicable in USco if the low disc fraction limit of ~20% is 
taken into account. In comparison, formation from very low- 
mass cores (1) would require higher turbulence and likely a 
denser region. Model (5) is also applicable in a dense region. 
A large fraction of the sub-stellar sources formed via these 
mechanisms could end up with discs, resulting in high disc 
fractions as found in the dense SOri or IC 348 regions. Thus 
while a dense environment could disrupt the discs around 
high-mass stars (e.g., Pfalzner et al. 2006), it could result in 
a high disc fraction for BDs. This can explain the large dif- 
ference in the BD and high-mass disc fractions in SOri and 
IC 348 (Fig, nil) . On the other hand, if mechanisms such as 
(4) are in play, the disc disruptions are minimal and all BDs 
formed should be able to retain their discs, irrespective of 
the stellar density in the region. This could explain the dif- 
ference in the BD disc fractions for Cha I and USco, which 
have similar stellar densities. 

The sharp decline in the BD disc fraction between 3 
and 5 Myr thus may not be an age effect but likely indicates 
that USco has formed its BDs differently than SOri or the 
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other young clusters. The one mechanism that seems the 
most apphcable in USco to explain the low BD disc fraction 
is the ejection model. All of the other models imply similar 
BD disc properties as found for the stellar members in the 
cluster. We do find the disc fraction for BDs in USco to be 
similar to the stars. Then perhaps any of the other formation 
mechanisms could work as well, although density could be 
an issue. We could argue that the disc fractions for the stars 
may also be indicative of the respective formation mecha- 
nism. A difference in BD and star formation mechanisms 
has been argued for in e.g., Thies & Kroupa (2007), and it 
has been suggested that stars and BDs are two correlated 
but disjoint populations with different dynamical histories. 
However, the high- and low-mass stars 'arrive' at this low 
USco fraction via a different track, i.e. via a transition to the 
debris phase. We have observational evidence of it in terms 
of the high debris disc fraction in USco for these higher mass 
bins, unlike the BDs. The similar disc fraction for stars and 
BDs in USco is probably not indicative of similar formation 
mechanisms for stars and BDs in USco, rather a coincidence 
that at this time stellar discs have evolved to a low fraction 
whilst BD discs started with a low fraction. The absence 
of any clear age dependency would explain the primordial 
nature of older BD discs. The ~10 Myr old TWA could be 
an exception where the extremely low-density environment 
might allow more BD discs to survive at early times of for- 
mation. A better argument could be that the disc fractions 
observed are a combination of evolution with age and the 
initial disc fraction due to the different formation mecha- 
nisms that are dominant in each cluster. 



6 SUMMARY 

We have studied the BD disc dissipation timescales in clus- 
ters at ages of ~1-10 Myr, and have investigated the various 
factors that may affect disc dissipation for BDs. Combining 
with several previous surveys, we find no clear age depen- 
dence for the BD disc fractions. The two most notable points 
are a sharp drop in the BD inner disc fraction between 3 and 
5 Myr, and the persistent primordial nature of these discs 
even at older ages. This indicates that while the overall BD 
inner disc fractions have decreased, the discs tend to remain 
in their primordial stages at ages as old as ~10 Myr. We do 
not find any convincing evidence that any of the well-known 
inner disc dispersal mechanisms are active in BD discs. We 
have instead considered the different BD formation theories, 
and have argued that USco may have formed its BDs differ- 
ently than SOri and other younger clusters. Thus the sharp 
decline in disc fractions may not be due to the difference 
in ages but may be indicative of the different BD formation 
mechanisms in these clusters, and there may not be a clear 
age dependence for the BD disc fractions. 

Our work also includes a SED classification scheme 
based on the Ks and WISE bands of 3.4-12fi.m. We find the 
WISE [3.4]- [12] color to provide a clear distinction between 
the photospheric and the disc population. 
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Figure 2. The WISE 2-12n.m slope plotted against spectral type for Taurus sources. Symbols are the same as in Fig.[T] The SpT of 0-9 
imply M0-M9 types. Dashed lines mark the thresholds used for SED classification. 
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Figure 3. ctyyiSE versus various WISE colors for Taurus sources. The [3.4]-[12] and _ft's-[12] colors show the strongest correlation v^ith 
the slope {top panel). Typical uncertainties in the colors are between 0.01 and 0.08 mag. 



Disc Frequencies for Brown Dwarfs 13 




3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 

Spectral Type 



Figure 4. The WISE 2-12|xm slope plotted against spectral type for our USco sample. The BD and M dwarf discs found in this work 
are indicated by red and blue symbols, respectively. The SpT of 3.5-8.5 imply M3.5-M8.5 types. 
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Figure 5. WISE and WISE-NIR ccds for the USco sources. The BD and M dwarf discs found in this work are indicated by red and blue 
symbols, respectively. The primordial discs from Carpenter et al. (2009) are indicated by green circles. Black symbols denote the debris 
discs from Carpenter et al. survey. The range in uncertainties are 0.03-0.06 mag for the Xs-[3.4] color, 0.03-0.15 mag for the [3.4]-[12] 
color, and 0.04-0.2 mag for the [4.6]-[22] color. In the bottom panel, only the sources with a detection in the 22|xm band have been 

plotted. 



Disc Frequencies for Brown Dwarfs 15 



^. 1.5 



CO 




G5 MO 

Spectral Type 




G5 MO 

Spectral Type 

Figure 6. Top: The WISE [3.4]-[12] color versus spectral type. Symbols are the same as in Fig. [5] Data for the earlier type discs is from 
Carpenter et al. (2009). Bottom The 12|j.m to Ks flux ratio versus the spectral type. We have also included the BD discs from Scholz et 
al. (2007) and Riaz et al. (2009). These are also denoted by red symbols. 
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Figure 7. The BD discs in USco. 
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Figure 8. The M dwarf discs in USco. 
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Figure 8. Contd. 
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Figure 9. Ha equivalent width versus the spectral type for the USco disc sources from the present work (black open squares). Also 
included are the BD discs reported in Scholz et al. (2007) and Riaz et al. (2009). Dashed line is the empirical accretor/non-accretor 
boundary from Barrado y Navascuos & Martin (2003). The SpT of 3-9 imply M3-M9 types. 
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5 

Age [Myr] 

Figure 10. Disc fractions versus stellar age for the 3 mass bins considered: top-high-mass, middle-low-mass, bottomSDs. Both Cha I 
and IC 348 have ages of ~2 Myr, we have plotted the IC 348 point at 2.5 Myr for clarification. Primordial disc fractions are denoted by 
filled symbols; open symbols for the high- and low-mass stars represent the debris disc fractions (Table [2]|. 
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Figure 11. Disc fractions versus stellar density for the Taurus (~1 Myr), Cha I (~2 Myr), IC 348 (~2 Myr), SOri (~3 Myr) and USco 
(~5 Myr) regions. Symbols and the figure sequence is the same as in Fig. 1101 For the high- and low-mass stars, the total disc fractions 
(primordial-l-debris) have been plotted. 
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Figure 13. BD disc fractions versus BD/Star number ratio. 
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44.5 


10.717 (0-022) 


49.9 


9.081 (0-037) 


29.7 


6.934 (0.083) 


13.1 


12.805 


(0.023) 


12.167 


(0.024) 


11.76 (0-021) 


-1 


31 




US00I6 


10 


05 


-19 


19 


36 


M6 


-49.8 


12.179 (0.024) 


44.3 


11.624 (0-024) 


44.5 


10.202 (0-067) 


16.3 


9.083 (0.531) 


2 


14.218 


(0.029) 


13.431 


(0.025) 


12.7 (0-033) 


-1 


42 




UScoie 


19 


40 


-21 


45 


35 


M6 


-36.5 


11.646 (0.026) 


42.2 


10.986 (0-021) 


50.9 


9.054 (0-034) 


31.7 


7.299 (0.11) 


9.9 


13.223 


(0.029) 


12.537 


(0.023) 


12.115 (0-024) 


-1 


09 




UScoie 


28 


47 


-24 


28 


14 


M6 


-182.7 


11.277 (0.024) 


44.7 


10.589 (0-023) 


47.5 


8.779 (0-055) 


19.7 


7.312 (0.197) 


5.5 


13.245 


(0.026) 


12.533 


(0.023) 


11.952 (0-023) 


-1 


06 




UScoie 


30 


27 


-23 


59 


09 


M6 


-34.7 


11.034 (0.026) 


42.1 


10.454 (0-022) 


48.6 


8.466 (0-046) 


23.5 


6.092 (0.074) 


14.6 


12.619 


(0.022) 


11.938 


(0.024) 


11.479 (0-021) 


-1 


12 




UScoie 


09 


59 


-23 


45 


19 


Me. 5 


-34.4 


11.364 (0.025) 


44 


10.977 (0-022) 


49.1 


9.994 (0-059) 


18.5 


8.41 (99.999) 


1.1 


12.635 


(0.024) 


11.996 


(0.024) 


11.624 (0-024) 


-1 


.8 




UScoie 


14 


33 


-22 


42 


13 


Me. 5 


-73.4 


12.904 (0.027) 


39.8 


12.286 (0-027) 


39.7 


11.405 (0-214) 


5.1 


8.511 (99.999) 


1.1 


14.232 


(0.03) 


13.664 


(0.033) 


13.24 (0-032) 


-1 


.7 




UScoie 


15 


14 


-23 


04 


26 


Me. 5 


-35.9 


13.487 (0.031) 


34.9 


12.826 (0-031) 


35.3 


11.412 (0-184) 


5.9 


8.701 (99.999) 


0.7 


14.814 


(0.042) 


14.193 


(0.032) 


13.915 (0-055) 


-1.40 




UScoie 


06 


04 


-20 


56 


45 


M7 


-159 


12.116 (0.027) 


40.2 


11.609 (0-025) 


43.5 


9.868 (0-058) 


18.8 


8.056 (0.211) 


5.1 


13.528 


(0.03) 


12.906 


(0.023) 


12.475 (0-027) 


-1 


34 




UScoie 


14 


21 


-23 


39 


15 


M7 


-43 


13.58 (0.033) 


33.3 


12.958 (0.032) 


34 


11.205 (0.147) 


7.4 


8.704 (0.387) 


2.8 


15.053 


(0.037) 


14.363 


(0.033) 


13.999 (0.056) 


-1 


24 




UScolS 


54 


20 


-21 


35 


43 


M8 


-29.4 


13.265 (0.029) 


37.5 


12.548 (0.029) 


37.8 


10.981 (0.126) 


8.6 


8.538 (0.316) 


3.4 


15.016 


(0.037) 


14.255 


(0.051) 


13.772 (0.049) 


-1 


25 




UScoie 


08 


30 


-23 


35 


11 


M8.5 


-97.4 


13.357 (0.029) 


37.8 


12.833 (0.031) 


34.9 


11.41 (0.19) 


5.7 


8.876 (99.999) 


0.4 


14.916 


(0.032) 


14.214 


(0.036) 


13.814 (0.054) 


-1.47 
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Table 2. Disc Fractions and Stellar Densities 



Name 


Age (Myri^ 


Distance {pc)^ 




Disc Fraction [%] 




Stellar Density [objs/deg^; objs/pc'^] 


BD/Star Ratio 








high-mass 


low-mass 


BDs 






Taurus'^ 


1-1.5 


140 


75±26 


60±20 


43±15 


28; 4.7 


'-'■^'-0.06 


Cha P 


2-3 


160-170 


65±22 


47±16 


55±19 


17; 2.2 


25+°-2* 


IC 348'= 


2-3 


315 


42±14 


43±15 


58±10 


415; 15.1 




'-'■^^-0.04 


SOri'^ 


3±2 




28±10 


38±13 


57±19 


345; 9.1 


0.23±0.02 


UScq£ 


5±2 


145±2 


4±3, 18±6 


19±6, 8±2 


26±9 


17; 2.6 


0.41±0.03 


TWi>3 




58.8±5.5 


50±20 


9±3, 18±6 


50±20 


0.06; 0.05 


0.16±0.04 



°' Reference for the median age estimates are: Taurus (Luhman 2004), Cha I (Luhman 2007), IC 348 (Luhman et al. 2003), SOri (Zapatero 
Osorio et al. 2002), USco (Preibisch & Zinnecker 1999), TWA (Barrado y Navascues 2006). 

^ References for distance estimates are: Taurus (Luhman 2004; Luhman et al. 2010), Cha I (Luhman 2007), IC 348 (Luhman et al. 2003), 
SOri (Ferryman et al. 1997), USco (de Bruijne et al. 1997), TWA (Biller & Close 2007). 

These discs are classified as primordial/evolved discs. The debris disc fraction for the younger clusters is negligible (~l-2%). The disc 
fractions for these clusters are from Luhman et al. (2008; 2010). 

The first value is for primordial discs, second value for debris discs. The BD discs are all classified as primordial/evolved discs. The 
disc fractions are from Carpenter et al. (2009), Riaz et al. (2009), Scholz ct al. (2007) and the present work. The total disc fraction for 
high-mass stars is 22±7%, and for low-mass stars is 27it9%. 

^ For the high-mass stars, all discs are classified as debris sources. For the low-mass stars, the first value is for primordial discs, second 
value for debris discs. The BDs are all classified as primordial/evolved discs. The disc fractions are from Low et al. (2005) and Riaz et 
al. (2009). The total disc fraction for low-mass stars is 21±7%. 
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APPENDIX A: COMPILATION OF DISC 
FRACTIONS AND STELLAR DENSITIES 

The disc fractions for Taurus, Chamaeleon I (Cha I), IC 
348 and a Orionis (SOri) are from Luhman et al. (2008) 
and Luhman et al. (2010). The age estimates for Cha I and 
IC 348 are ~2 Myr (e.g., Luhman et al. 2008), and 3±2 
Myr for SOri (e.g., Zapatero Osorio et al. 2002). Luhman et 
al. (2008; 2010) have converted the spectral types to stel- 
lar masses by combining the temperature scale of Luhman 
et al. (2003a), and the evolutionary models from Baraffe et 
al. (2003). They have considered the mass limits of 0.1 Mq 
and 1 Mq at SpT of M6 and K5, respectively. In Fig. \Wi 
the sub-stellar disc fraction is for sources with SpT between 
M6 and M9, implying masses between ~0.1 and 0.02Mq, 
the low-mass disc fraction is for SpT of K5-M5 (O.I-IMq), 
and the high-mass disc fraction is for SpT earlier than K5. 
The earliest type in these clusters is a mid-B star (B5 or 
B6; Luhman et al. 2003a; 2007). From the models of Seiss 
et al. (2000), at an age of 2-3 Myr, a B5 SpT corresponds 
to a mass of ^4 Mq. The high-mass disc fractions plotted 
are thus approximately for the mass range of 1-4Mq. Since 
most of the SOri sources do not have spectral classification, 
Luhman et al. (2008) have used the Mj magnitudes as a 
proxy for stellar masses, and have considered the mass lim- 
its of IMq and O.IMq at Mj of 3 and 7 mag, respectively. 
From the Mj ranges listed in Luhman et al. (2008; Table 4), 
we have considered Mj between and 3 mag for high-mass 
stars, Mj=4-8 mag for low-mass stars, and A4"j=8-10 mag 
for sub-stellar sources in SOri. These Mj ranges correspond 
to similar mass ranges as considered for the other three clus- 
ters. In USco, the BD disc fraction (26±9%) is the combined 
disc fraction from the work of Scholz et al. (2007), Riaz et 
al. (2009), and the present work. The disc fraction for low- 
mass stars in USco is for the K5-M5 stars from Carpenter 
et al. (2009), and includes the 36 M dwarf discs (M3.5-M5) 
found in the present survey. The disc fraction plotted for 
high-mass stars in USco is for SpT earlier than K5 obtained 
from Carpenter et al. (2009). In the ~10 Myr old TWA, the 
BD disc fraction (M6-M8) is from Riaz et al. (2009), and the 
disc fraction for the low-mass stars (K5-M5) is from Low et 
al. (2005). Among the higher mass stars in TWA, the Low 
et al. sample consists of a G5 star (TWA 19A) and an AO 
star (HR 4796A). TWA 19A is found to be photospheric, 
while HR4796A shows excess emission in the Spitzer/MlPS 
bands. This results in a 50±20% disc fraction for the high 
mass stars. We note that the TWA is a sparsely populated 
association with a small number of 8 discs found to date. 
The disc fractions for this association thus have large un- 
certainties. The error bars for the disc fractions are the l-cr 
Gaussian distribution uncertainties. 

The quoted disc fractions in other clusters are based on 
deep Spitzer surveys, which have high sensitivities capable 
of detecting the photospheres for the brown dwarfs. If we 
consider the Spitzer 8|i.m band, which is the closest in wave- 
length to the WISE 12^j.m band, then the predicted photo- 
spheric flux at this wavelength for an M6 dwarf is ~0.3 mjy, 
and for an M8 dwarf the flux is ~0.13 mJy. The deep sur- 
veys conducted with Spitzer in SOri, Taurus, Cha I and IC 
348 have a l-cr sensitivity of ~0.01-0.02mJy at 8(i.m. These 
surveys were thus capable of detecting the photospheres and 
therefore the smallest excesses possible. The smallest 12(j.m 



excess that we can detect with WISE is Fobs/Fphot ~ 3-5 
for M6-M8 objects. We have checked in the WISE catalog 
for detection of all of the Spitzer known discs, in order to de- 
termine the fraction of discs that WISE might have missed 
around M6-M8.5 type objects due to limited sensitivity. In 
SOri, all of the Spitzer brown dwarf discs at masses of ~0.08- 
0.04 Mq, corresponding to spectral type of ~M6-M8, were 
recovered at S/N > 3 in WISE 1, 2, 3 bands. The S/N is 
~l-2 in the 22vim band. For lower masses of ~0. 04-0. 02 Mq 
(M8-M9), none of the sources were detected in WISE 12(i.m 
or 22nm bands (S/N < 1), and detections were only obtained 
in the 1 and 2 bands at S/N~10. From the number of missed 
discs, we estimate that the WISE BD disc fraction in SOri 
would be --^10% lower than the Spitzer fraction. SOri is also 
the farthest cluster (~352 pc) considered in our study, which 
explains why the later type discs were not detected. IC 348 
is another distant cluster (~315pc) and WISE has missed 
the M8.25-M8.75 discs in this cluster, which results in a 
BD disc fraction which is -^8% lower than that estimated 
by the Spitzer survey. We note that due to the 3-a quality 
criteria, the full sample detected in the WISE 1, 2, and 3 
bands at S/N^3 would actually be smaller in size than the 
original sample, resulting in a higher disc fraction. Combin- 
ing these limitations, we estimate that the WISE-recovered 
disc fractions in SOri and IC 348 will be ~52% and ^53%, 
respectively, which are similar to the Spitzer fractions of 
58% and 57% (Table [IJ. In the Taurus and Cha I regions 
(~140-160pc), all of the brown dwarf discs as well as the 
full samples down to a spectral type of M9 were recovered 
by WISE, at a ^3 sigma level at 12(j.m, as well as in WISE 
bands 1 and 2. Therefore the WISE-recovered disc fractions 
in these regions are similar to the Spitzer ones. The WISE 
data on TWA (--^GOpc) has not been released yet, but from 
the Spitzer/IRS data discussed in Riaz et al. (2008), the ob- 
served flux density at 12(i.m is ~5 mJy for the TWA brown 
dwarf discs, which is higher than the 5-a sensitivity at this 
wavelength (~0.7 mJy). These discs thus should be recov- 
ered with WISE. We therefore do not expect any significant 
changes in the overall trend with age observed for the BD 
disc fractions shown in Figs. 1101 and the USco disc fraction 
would still be a factor of ~2 lower then the SOri fraction. In 
Figs. 1101 [TT]and Table [51 we have used the Spitzer fractions 
quoted in the literature, and not the WISE-recovered disc 
fractions. 

It is important to determine the wavelengths at which 
these disc sources show excess emission. The disc fractions in 
Taurus, Cha I, IC 348 and SOri have been determined based 
on observations made in the Spitzer/IRAC bands of 3.6- 
8nm and MIPS 24(j.m band. The intermediate wavelengths 
of ~8-24[j.m are important in making a distinction between a 
primordial and a debris disc. Luhman et al. (2010) have clas- 
sifled sources that show excess emission at 24(j.m but photo- 
spheric emission at shorter wavelengths as debris discs, while 
primordial discs are sources that show excess emission at 5.8, 
8 and 24(i.m. Intermediate between the primordial and debris 
discs are the 'evolved' discs that show weaker excess emis- 
sion in these bands compared to the primordial discs. These 
evolved discs are similar to the primordial transition discs 
discussed in Section i]4.1l Based on this definition, there are 
just 4 debris discs in Taurus, 2 in IC 348, 1 in Cha I, and 
2 in SOri (Luhman et al. 2010). We can thus consider that 
the fraction of such debris systems with excesses at ^24(i.m 
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are negligible in these clusters (~l-2%). A small ~12% of 
the discs in Taurus are classified as evolved discs, while a 
large fraction is found to be primordial. Similar is the case 
in the other young clusters of IC 348, Cha I and SOri, where 
nearly 85% of the disc population is found to be primordial, 
and ~10-15% are evolved sources (Luhman et al. 2010). 

In USco, 19 of the 23 earlier type discs (SpT<K5) show 
excess emission at 24fi.m only, and have been classified as de- 
bris discs by Carpenter et al. (2009). The rest of the 4 discs 
have excesses at 8 and 24(j.m and have been classified as pri- 
mordial discs. Among the K5-M5 group, 58 are primordial 
discs with both 8 and 24fi.m excesses, while 10 are debris 
with 24(im excess but weak or no detectable 8[im excesses. 
The Carpenter et al. criteria is thus consistent with Luh- 
man et al. (2010), wherein a debris disc is classified based 
on 24)i.m excess only, while a primordial disc will show ex- 
cess emission at both 8 and 24(j.m. Carpenter et al. have not 
considered the intermediate evolved disc category, but the 
primordial discs in their survey are found to show weak or no 
excess shortward of 8\Lm, so these are likely to be evolved or 
primordial transition discs. In our present USco work based 
on the WISE bands, it is only the 12[i.m point that can make 
such a distinction. All of the discs from this survey have an 
excess at 12jim, and we do not have any discs that show 
an excess at 22jxm only. Therefore none of these have been 
classified as debris discs. As noted in Section tj4.1l the col- 
ors for these discs are also similar to the primordial discs 
classified by Carpenter et al. (2009). Only six of these show 
small excess at 3.4 and 4.6(i.m, so a majority are primordial 
transition discs. For the other BD surveys in USco, all of the 
discs from Scholz et al. (2007) show excess emission at 8vim 
and at longer wavelengths. The discs from Riaz et al. (2009) 
also have excess at 8|J.m and at longer wavelengths, with one 
disc showing excess in the IRAC shorter wavelengths bands 
of <8(i.m. Thus all of the USco BD discs would fall into the 
primordial/primordial transition disc category. 

In the TWA, the BD discs have excesses at ~5-24|im, 
and can be classified as primordial (Riaz & Gizis 2008). 
These discs show weak signs of inner disc clearing. Among 
the low-mass stars, TW Hya and Hen 3-600 can be classi- 
fied as primordial transition discs, since both show strong 
excesses at 10^-m and longer wavelengths, but photospheric 
emission at shorter wavelengths. The sources TWA 7 and 
TWA 13 show excesses only at 70 (im, while two other discs 
show weak or no excess shortward of 24|i.m. These 4 discs can 
be classified as debris systems. In Fig. 1101 we have plotted 
the total (primordial/evolved -|- debris) disc fraction. It is 
difficult to clearly separate the primordial and evolved disc 
sources since the definitions vary from survey to survey. As 
noted earlier, the debris disc fraction for the younger clus- 
ters is negligible (~l-2%), therefore a majority of these can 
be considered as primordial discs. In comparison, a major- 
ity of the high- and low-mass discs in USco and TWA are 
in the debris phase. In Table [2l we have listed the separate 
primordial and debris disc fractions for these clusters. 

Figure [11] plots the disc fractions versus the stellar den- 
sity in a cluster. By stellar density, we imply the total num- 
ber of objects (stars -I- BDs) detected in a given survey area. 
The stellar density in Taurus is from the work of Luhman 
(2004), where 112 objects (96 stars and 16 BDs) were con- 
firmed as members from imaging of a 4 deg^ region, im- 
plying a stellar density of 28 objects per square degree. In 



Cha I, Luhman (2007) detected 85 objects (70 stars and 15 
BDs) in his IK2 survey covering a 1.25° radius (area=4.9 
deg^), implying a stellar density of ~17 objects per square 
degree. In IC 348, Luhman (1999) found 288 members (in- 
cluding 23 BDs) in a 50'x50' central region of the cluster, 
implying a stellar density ~415 objects per square degree. 
In SOri, Lodieu et al. (2009) found 271 objects (198 stars 
and 73 BDs) in a 30' radius from the cluster center, im- 
plying a stellar density of ~345 objects per square degree. 
In USco, Lodieu et al. (2011) found 81 stars and 31 BDs 
in a 6.5 deg^ region surveyed, resulting in a stellar density 
of ^17 objects per square degree. The TWA, in compari- 
son, is a loose association of about 30 confirmed members. 
From Webb et al. (1999), the TWA searches were mainly 
conducted in a 12° radius of the approximate cluster center. 
This implies a stellar density of ~0.06 objects per square 
degree. We note that there are several other surveys con- 
ducted for these clusters/associations, which would result 
in slightly different stellar densities. For example, in Tau- 
rus, Luhman et al. (2003b) found 92 members in a 8.4 deg^ 
area, implying a stellar density of ~11 per square degree. 
In USco, Slesnick et al. (2008) detected 1000 objects in a 
62 deg^ area, implying a density of 16 per square degree. 
In IC 348, Luhman et al. (2003a) found 122 members (in- 
cluding 15 BDs) in a 42'x28' region, resulting in a density 
of ~373 sources per square degree. Our results for stellar 
densities indicate SOri and Cha I to be denser than Taurus 
and USco, which is in agreement with what has been pre- 
viously noted by several surveys. In Tabled we have listed 
the stellar densities plotted in Fig. [TT] The stellar densities 
are listed in both stars/deg^ and stars/pc^ units. 

In Fig. 1131 we have plotted the disc fractions versus 
the ratio of number of BDs to the number of stars in these 
clusters. The BD/Star number ratios are listed in Tabled 
The ratios for Taurus, Cha I and IC 348 have been taken 
from Andersen et al. (2008), the ratios for SOri and USco 
are from Lodieu et al. (2009; 2011), and the ratio for TWA 
has been calculated using data from Low et al. (2005) and 
Riaz & Gizis (2008). Note that the mass range in Andersen 
et al. (2008) are different than Lodieu et al. (2009; 2011). 
Andersen et al. have used 0.08-1.0 A4q for "stars", while 
the range is 0.08-0.35 Mq in Lodieu et al. (2009; 2011). For 
BDs, the mass range is 0.03 - 0.08 Mq in Andersen et al., 
and 0.02 - 0.08 Mq in Lodieu et al. However, the ratios are 
consistent within the uncertainties from these surveys. The 
values provided in Andersen et al. (2008) have been calcu- 
lated using data from Luhman et al. (2003) and Luhman 
(2004; 2007) surveys. The mass ranges are slightly different 
but the ratios are consistent within the uncertainties. The 
IC348 sample from Luhman et al. (2003) is unbiased down 
to 0.03 Mq and not 0.02 Mq as used by Andersen et al.. 
However, BD/star ratios in Luhman et al. (2003) work are 
defined down to 0.02 Mq, which is also the limit Lodieu 
et al. (2011) have used. The Taurus sample form Luhman 
(2004) work is unbiased down to 0.02 Mq. In IC 348, Ander- 
sen et al. have listed a ratio of 0.12 (0.08-0.18), Luhman et 
al. (2003) give a ratio 0.18±0.04, and Lodieu et al. (2011) list 
a ratio of 0.20±0.03. These ratios are all consistent within 
the uncertainties. 

Looking at other BD/Star ratio values for the regions 
listed in Tabled we have found the following: 
1) Taurus: 0.17 (0.11-0.25) from Andersen et al. (2008) and 
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0.14±0.04 from Luhman (2004) 

2) SOri: 0.23±0.02 from Lodieu et al. (2009) 

3) Cha I: 0.25 (0.13-0.53) from Andersen et al. (2008) and 
0.26±0.06 from Luhman (2007) 

4) IC348: 0.12 (0.08-0.18) from Andersen et al. (2008) and 
0.12±0.04 from Luhman et al. (2003) 

5) USco: 0.41±0.03 from Lodieu et al. (2011) 

In Table [5] and Fig. 1131 we have used the values from 
Andersen et al. (2008) for Taurus, Cha I and IC 348, and 
from Lodieu et al. (2009; 2011) for SOri and USco. 



